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SYNOPSIS 


Blexural and torsional natural frequencies of a sh aft , 
which is loaded and supported periodically, are, investigated 
in this thesis. The analysis is carried out by using a ’wave 
approach*. In this approach, first the free harmonic wave 
propagation, hoth in torsion and in flexure, in an otherwise 
identical infinite shaft is studied. The results are then 
utilized to compute the natural frequencies of the finite 
shaft. Phis method renders the computational effort indepen- 
dent of the number of spans in the shaft. 

By proper choice of parameters, it is shown that 
operating speed ranges, free from both bending and torsional 
resonances, can be obtained. Suggestions for widening such 
speed ranges of operation are also included. 



CHAPTER 1 


INTRODUCTION 


1.1 Introduction 

Most of tlie marine vessels use steam and gas turbines 
for both propulsion and power generation. In these, the 
shafts for transmitting power need a safe speed range of 
operation without any resonance. It is also important to 
know the vibration characteristics of a ship’s diesel 
installation consisting of a stepped shaft on multiple 
supports with a series of discrete inertias. Furthermore, 
there are instances where the shaft rigidity is of primary 
consideration as compared to its strength, e.g., the cam- 
shaft in various engines. These shafts are also provided 
with a large number of supports in order to minimise the 
shaft deflection. This is essential for proper functioning 
of the cam-follower system. 

All these shafts may be excited by loads which 
generate vibration both in torsional and flexural modes. 

These vibrations, if excessive, may give rise to mal- 
functioning and even fatigue failure in some cases. A shaft with 
too many lumped inertias has closely clustered torsional 
natural frequencies, likewise, a shaft on a large number 
of supports has closely clustered flexural naimral frequen- 
cies. The number of frequencies in every cluster depends 
on the number of lumped inertias and the supports. 
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respectively. Moreover, the lending and the torsional natural 
frequencies may also overlap. All these result in a very 
narrow range of operating speeds for such shafts. Hence, 
while designing these shafts, it is essential to have an 
accurate knowledge of all the natural frequencies. 

To this end, as a first choice to a designer, the 
shaft may he idealised as one loaded and supported periodi- 
cally. After determining the natural frequencies the study 
can he oriented towards an attempt to put the bending and 
the torsional natural frequencies in a narrow range. This, 
if achieved, may result in a wider range for the operating 
speed. 

1.2 Short Review of Previous Work 

Extensive work is reported on finding the torsional 
frequencies of crank shafts of multi-cylinder engines [1,2] 
where crank shafts are idealised as lumped rotors connected 
hy inertialess shafts. Such a model is accurate as long as 
the shaft inertia is reasonably smaller than the smallest 
inertia of the rotors [ 3 ]. Natural frequencies of a heavy 
shaft, e.g., that of a heavy marine engine with rotors, were 
developed hy Lewis [4]. He developed a graphical technique 
of finding out the frequencies of a single— mass elastic 
system and later, extended his method to systems of several 
degrees of freedom. The frequency equation for the special 
case of a single shaft with rotors at either end is available 
in Ref. [5], while for the general case, a hand-computational 
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method is outlined by her Wilson [l], In reference [l], 
the torque summation for a system performing free torsional 
vibration put equal to zero for finding out the exact 
natural frequencies. The difficulties often encountered with 

t 

the torque summation method (Holzer’s method) are (l) in 
estimating the initial value of Ob 2 , and (2) in selecting 
the second trial value ofg 2 if the initial trial value 
does not satisfy the torque summation equation. The 
estimation of the initial trial value depends on judgement 
and experience, hut reasonable trial values for the lower 
modes can be obtained by reducing the multi-disc system to 
a two— or three-disc systems. Grouping together discs that 
have shafts with relatively high stiffness between them, or 
neglecting discs with inertias that are small compared to 
the other discs, however, would give a fair approximation, 

Maltbaek [6,7] studied free torsional vibrations of a 
stepped shaft with a series of discrete inertias and 
torsional springs acting on the shaft. In this analysis, 
a determinant of order (2nx2n) is solved to obtain the 
natural frequencies, where n is the number of discrete 
inertias. He also described [6] in detail the influence of 
a discrete inertia (its value and position) on the free 
torsional vibration of a uniform shaft with torsionally 
elastic end constraints. Rao [8] derived a simple n'*'* 1 
order frequency determinant to compute the frequencies and 
modes of an n stepped shaft with rotors. As an example, 
the method was used to find the frequencies and modes of a 
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heavy homogeneous engine in which the effect of shaft inertia 
was shown to be small for higher modes. Natural frequencies 
of such systems decrease linearly with the shaft inertia 
parameter when the shaft is very heavy. 

It is seen from the above literature that the analysis 
of long shafts with too many rotors becomes cumbersome and 
unwieldy by the conventional methods. If the shaft is 
loaded periodically , a’wave approach' can, however, be used 
profitably. In this method, the amount of computation 
remains independent of the number of lumped inertias. In 
tlao case of flexural vibration of periodically- supported 
beams, wave propagation theory has been applied to find the 
modes and natural frequencies [9-14]. 

The concept of wave propagation in periodically- 
supported, undamped structures was introduced by Heckl 
[9,10]. Later Mead [ll] analysed the free harmonic wave 
propagation in an infinite beam supported periodically by 
both rigid and flexible supports. He obtained distinct atte- 
nuation and propagation zones of free wave motion in such a 
beam. Using this wave approach, Sengupta [12] reported a 
graphical method of finding natural frequencies of periodi- 
cally-supported finite beams with different end conditions, 

1,3 Objective and Scope of the Present Work 

The major objective of the preset t work is to study 
harmonic torsional wave propagation in an infinite shaft with 
periodic rotors, taking into account of nhe inertia of the 



shaft, file rotors are idealised as circular discs. The 
discs and shafts are assumed to he of the same material. 

The complex propagation constant for such a shaft is 
obtained as a function of frequency. Alternate propaga- 
tion and attenuation bands of the free wave motion is 
plotted against the frequency, from these plots, natural 
frequencies of finite shafts are obtained. 

Next, to consider the bending vibration of such a 
periodically loaded shaft, it is assumed that the shaft 
is supported by identical, equi-distant bearings, further- 
more, only one. rotor is supposed to ha located between two 
consecutive bearings. Short bearirg s are idealised as line 
supports whereas long b .arings are considered to provide 
rotational constraints to the bending vibration. Initially, 
all the bearings are considered to be rigid so far as 
transverse oscillation is considered. The bending natural 
frequencies are also obtained by using the wave, approach 

[ 13 ]. 

In the third phase, an attempt is made to match the 
groups of bending and torsional frequencies as far as 
possible. To this end, a parametric investigation is 
carried out by varying the parameters L/d, D/d, ^ and k r j 
where L is the bearing spacing, D and d refer respectively, 
to the disc and shaft diameters, ^ represents non-dimensional 
distance of the disc from the adjacent bearing on the left 
and £ r indicates the non-dimensional rotational stiffness 
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of the bearing. With a good matching of these clusters of 
bending and torsional frequencies, a wider operating speed 
range for the shaft could be available. In such a range the 
shaft will be free from both bending and torsional oscillations. 
At the end, a short discussion is also added wherein 
the bearings are considered to be flexible in the transverse 
direction, ho detailed analysis with flexible bearings is 
included in the present work. She practical limitation of 
such an idealised shaft which is both loaded and supported 
periodically should be kept in mind. However, the results 
obtained from the analysis presented in the thesis can act 
as a good guide to a designer for his first choice. 


CHAPTER 2 


FREE TORSIONAL AND BENDING VIBRATION OF A SHAFT LOADED AND 

SUPPORTED PERI ODI CADDY 

2.1 Theory of Harmonic Wave Propagation 

2.1.1 The concept of wave propagation 

In an infinite shaft or beam tjnpe periodic structures, 
harmonic waves can propagate freely (without decaying) in 
certain frequency ranges called propagation zones. Frequen- 
cies outside these propagation zones are called attenuation 
zones. In these zones, a wave once started will de-cay as 
it spreads outwards. The complex propagation con st ant 41 which 
characterises the wave motion is the measure of the rate of 
decay and the change of phase of the wave motion over the 
distance between adjacent supports. There always exist 
positive and negative propagation constants (+ p). The 
negative sign in front of the complex number implies that the 
wave is propagating, or decaying, from left to right and the 
positive sign implies the reverse direction. The propagation 
constant is a function of frequency and system parameters. 

2.1.2 Evaluation of propagation constants 

A periodic system consists of a number of identical 
elements, E, each having a single coupling co-ordinate with, 
respect to its neighbour at either ends (Fig. 2.1(a)). Such 
a system is referred to as a mono-coupled system. Each 
element E is characterised by its receptances, which relate 
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(b) 


FIG- 2-1 (a) BLOCK DIAGRAM OF A SYSTEM OF 
MONOCOUPLED ELEMENTS 

FIG-2-1 (b) FORCES ON, AND DISPLACEMENTS 
OF A SINGLE ELEMENT ' 
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the harmonic displacements at the coupling co-ordinates to 
the harmonic forces at those co-ordinates, let the harmonic 
displacements at the left and right hand coupling co-ordinates 
of a single element he q A exp(iOOtt) and q B exp(i& 3 t) and the 
harmonic forces acting opposite to one another at these co- 
ordinates he P A exp(icat) and F B exp(ic^t), respectively 
(Fig. 2 . 1 (h)). The displacements and forces are related as 
follows through the receptance matrix [8 ] s 




( 2 . 1 ) 


where is the harmonic displacement at co-ordinate i due 
to a unit harmonic force at co-ordinate 3 . The receptances 
Pjj’s are functionsof frequency and the element characteristics. 

If an infinite number of elements like ! E* are coupled 
together, free motion of the system is possible when the 
followiig recurrence relationship holds good for every element: 


q B = q A 


( 2 . 2 ) 



e 


(2.3) 


For a linear system P AB = - P BA , so that there are three 
significant receptances. Further, for a symmetric element 

fhpsp r'oflur'.p -hn -twn as B . , = — 
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Using Eqns. 2,2 and 2,3 , the pair of equations 2.1 
hecome 


lA = «AA + M *A 
and 4 A = (P BA + e^ 1 P BB ) Jq 

for non- trivial solution of, E^ and Eg, in. case of a 
system; > one gets 

“ ^BB + (e^ + e ^3 P AB = 0 


hence 


C0S11 [1 = 


^BB ~ $ m ) 


2P 


AB 


Eqn, 2.6 simplifies to a quadratic equation in e 1 ” 1 as 


e 2 ^ + C(P M - PbbVPabI ^ + 1 = 0 




the solution of which is 


(2,4) 


(2.5) 

linear 


( 2 . 6 ) 

(2.7) 


( 2 . 8 ) 


e 


^1,2 


r + 



(2.9) 


where r = (p BB - = cosh V 


( 2 . 10 ) 


2.1.3 Dependency of propagation constant on frequency 

Considering an undamped periodic 'system, the term r in 
Eqn, 2.10 is entirely real and varies between +cxl as the 
frequency is varied. If it is greater than + 1, p is 
entirely real (p r ) and there is a constant exponential rate 
of decay of the wave motion from one element to the next. 
Since p has no imaginary part, there is no phase difference 



between the motion in adjacent elements and the wave is of 
non-propagating form. If r is less than - 1, then p has 
the form p r +i7i;, the imaginary part implies that adjacent 
elements vibrate exactly in counter-phase and the wave is 
still of non-propagating type. These frequency ranges, for 
which p r ^ 0, are called attenuation zones. On the other 
hand, if r is real and is between -1 and +1, p is entirely - 
imaginary (ip^) which implies non-zero phase difference 
between the motions in adjacent elements, so that the wave 
is now propagating. There is no decay of the wave motion 
between successive elements as the real part of p is zero- 
The frequency ranges for which p r = 0 are referred to as 
propagation zones, 

2,1,4 Natural frequencies of finite periodic systems 

Wave propagation can occur provided p is imaginary, an 
this occurs when -1 < cosh p < +1. Thus the bounds of the 

r Ni 

frequency in which wave propagation takes place are gov erne 1 
by 

coshp = “ ^AA^^^AJB = i (2,11) 

T'or a symmetric element, (3^ = -p-gg, so that cosh p = 
“^Af^AB* Bhe bounding frequencies of the propagation bojilnds 
in this case are therefore that at which = + PaB* 

let us consider a single system element, disconnected 
from its neighbours, so that its ends are free. At its 
natural frequencies the receptances have infinite values, 



"but the ratios and Pjy^/P^g remain finite and determinate 

In particular, for a symmetric element they become -1 and +1, 
respectively, the sign of the later depending on the nature 
of the system and the mode. Hence at the natural frequencies 
of the free symmetric element, the finite receptance ratios 
make cosh p for the periodic system equal to + 1. 

If the same single element is locked at both of its 
ends, it has different natural frequencies. These are obtained 
from the pair of following equations : 


^A 

" ^AA P A 

+ ^AB E B 

2.12(a) 

^B 

= ^BA P A 

+ ^BB ^B 

2.12(b) 


lor a locked element, q^ = q B = 0 and the natural frequencies 
are the roots of the determinant equation 


P AA P AB 


P BA P BB 


0 


(2.13) 


As mentioned earlier, for a symmetric element putting 
- , and = - Pgg equation 2,13; reduces to 

"dj-eL ■ 0 

or ^aa/^ab = ± 1 


p 


BA “ 


(2.14) 


At the frequencies which satisfy Eqn, 2.14 (the natural 
frequencies), cosh p is again found to to + 1. It follows, 
then, that the bounding frequencies of a periodic system of 



symmetric elements (at which, cosh p = + 1) are also natural 
irequencies of the single element with either free ends or 
locked ends. 

If the element is unsymmetric, ^ However } 

the natural frequencies of the locked element are still the 
roots of Eqn, 2.13 and this may "be expressed in the form 

Pju/PaB = “ (2.15) 

Hence, at these frequencies 

cosh p - (PbB ” 

= \ ^“^AA^AB^ “ ^AB^AA^ 

But cosh p = \ (e^ + e"^), and so e^ = - f^j/^AB or ^ Li: 

These "both have real values and for the general unsymmetric 
element, are not equal to + 1 at the natural frequencies. 

Hence, p must have a non-zero real part and the natural 
frequencies of the locked, unsymmetric element must lie in the 
attenuation zones. 

In case of a single unsymmetric element with free ends, 
it can "be shown [ 13 ] that at its natural frequencies, {3^ = 

p 

- ^£ 5 * Once again, then, at the natural frequencies P^a/P *33 = 

- P_£b^BB 331 ^ so " i: - i:iese frequencies must lie in attenuation 
zones. Hence the natural frequencies of both free and locked 
unsymmetric elements lie in attenuation zones and do not coincide 
with the "bounding frequencies. 
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Considering a system of positive-definite aLements, it 
can "be proved that as the frequency approaches zero, propa- 
gation cannot occur, i.e. the lowest frequency range is an 
attenuation zone. Propagation begins at (or just above) the 
lowest natural frequency of the positive— definite element, 
which, is free at its ends. If the ’’free' 1 and "locked” 
natural frequencies of the single symmetric element occur 
alternately (as for certain uniform beam elements) the 
free natural frequencies constitute the lower frequency 
bounds to the propagation zones, whereas the locked frequen- 
cies constitute the upper freqiiency bounds [l 3 ]. 

If the system is semi-definite, propagation does 
occur as the frequency approaches zero, i.e. the lowest 
frequency range is a propagation zone. Propagation ceases 
at (or just below) the lowest of the natural frequencies 
of the semi-definite element with its ends either free or locked. 
This frequency is therefore an upper bound frequency. 

Now we consider a finite mono-coupled periodic system 
of N elements being excited by the force Pp at the left end 
P (Pig. 2.2). Since the system is mono-coupled, the total 
motion generated can be represented by a single positive- 
going characteristic wave and a single negative-going wave. 

If the contributions to the total displacement at P from 
these two waves are qp+ and qp- respectively, then 

q.p = q.p+ + Ip- (2.16) 
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FIG-2-2 BLOCK DIAGRAM OF A FINITE PERIODIC SYSTEM 
OF N ELEMENTS AND ARBITRARY BOUNDARIES 



FIG- 2-3 THE MULTIVALUED PROPAGATION CONSTANT 
CURVES VERSUS FREQUENCY . 
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At the boundary Q the displacement of the positive wave is 
4q+ = V 

and that for the negative wave is 



Thus, the total displacement at Q is 

1q = Hq+ Hq— = lp+ ® ^ + ip— ® ^ (2.I7) 

Defining the characteristic wave receptance (P ) as the 
harmonic response at one of the coupling co-ordinates per 
unit internal harmonic force at that co-ordinate one gets 
from Eqn. 2.4 

P W = = ^AA + ^ ^AB ( 2 , 18 ) 


Hence ,the characteristic receptances of positive and negative 
going waves are, respectively, given by 


fV = f>AA + Pjffl 2 - 19Ca) 

and = P^ + P^g 2. 19^) 

The total force at each end of the system shown in Dig. 2,2 
is given by 

Ip = ,1 P +/P W + + lp-/P w - 2.20(a) 

and Eq = lp+ e + Ip- e ^/^w” 2.20(b) 
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Supposing that the system is free at end Q such that Pq = 0,, 
equation 2.20(b) then yields 

9p- = ~(P W -/P W +) ? P + e - ® 11 (2.21) 

On substituting this, and Eqns. 2.19(a) and 2.19(b), into 
Eqn, 2.20(a) the direct receptance at end P is found to be 

Ppu = IpAp = P AA - feinh(U-l)p/sinh Up] 

.... 2.22(a) 

or Ppp = $AA ~ ^AB - cosl1 P “ cosh Up sinh p/sinh Up] 

.... 2.22(b) 

The transfer receptance, pQp = - PpQ, is likewise found to 
be 

PqP = Iq/^p = iVsinh Up] (2.23) 

These two receptances are entirely real, whether p is 
real or imaginary, i.e. whether the frequency is in a 
propagation zone or in an attenuation zone. The natural 
frequencies of the system which is free at both ends P and 
Q are those frequencies at which Pp P — >oc » These occur 
(a) when and ("M when sinh p/sinh Up — . 

Condition (a) occurs in general, in an attenuation zone 
but will occur at the bound of a propagation zone if the 
element is, symmetric. Condition (b) can only occur in 
a propagation zone, for then p is purely imaginary (ipp) 
and sinh p/sinh Up = sin pp/sin Upp. This becomes 
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infinite, whenever p^M = mm, provided ^ ^ 0 or + i. It 
has "been shown [l2] that the integer m is not equal to zero nor 
a multiple of 1, i.e. l^m<l, giving (1T—1 ) discrete values of 
M-jL 311 ^ correspondingly (IT— 1 ) discrete natural frequencies in 
the propagation hand. 

figure 2.3 shows the imaginary part of p (p^) plotted 
against frequency for the first propagation hand of a multi-- 
supported uniform beam. p^ is a multi-valued function of 
frequency since p ± = cos -1 [ (P BB - ** follows 

that each section of the cir ve between p^ = 0 and n is a 
mirror image of the section above it. 

let the finite system has four elements. It will 
resonate at the frequencies given by p^ = mn/4 (m ^ 0 or 
a multiple of 4) and these are the frequencies of the points 
marked a, b, c in fig. 2.3 for m = 1, 2, 3 and the points 
marked e, f, g for m = 5, 6, 7 etc. Since BC is the mirror 
image- of BA, the frequencies of the points a and g, b and f, 
and c and e, respectively, are identical. Thus it can be 
concluded that a finite periodic system with free ends and 
consisting of N mono-coupled unsymmetrie, positive-definite 
elements has 1-1 natural frequencies in each propagation band, 
corresponding to the frequencies at which p^ = mn/l (m = 1 
to 1-1), Close to, but just below each propagation band there 
is another natural frequency, corresponding to ft & A 

Mien the finite system with 1 elements is fixed at the 
right-hand end such that 1 q = 0, the direct receptance at 
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the 1 ef t-hand end is given "by 


p 1 ^AB ~ PaA 

Ppp - 

- PjLB (tT+l) p/sinh Np 


(2.24) 


Tile natural frequencies of a system oscillating in fixed-fixed 
mode (for which Ppp =0) are those at which either (a) 

(P AB “ ^AA ^BB^ = ^ or ? 0 °) sinh Up = 0, provided sinh(N+l)p ^ 0. 
Condition (a) is the equation for the natural frequencies of 
the single element, fixed at "both ends (Eqn. 2.13) * Bor a 
system, of general unsymmetric elements these frequencies lie 
outside the propagation "bands. If the elements are symmetric, 
they coincide with the upper "boundary frequencies of the 
propagation zones. Condition (h) can only "be satisfied in a 
progapation "band at those frequencies for which sinh Np = 
sin Np^ = 0 and sin (N+l) p^ ^ o, i.e. p^ = mrt/N, m = 1 to 
(N-l), Hence, these are the same (N-l) natural frequencies 
already discussed for the ’free-free 7 system, within each 
propagation zone. 


2.2 Evaluation of Receptances for a Periodic Shaft ELement 

Figure 2.4 shows a periodic shaft element of length 1 
with a disc of moment of inertia Ip at ^ ^ , where £ = ^ with 

x representing the distance of the disc from end A. 

The equation of motion for the shaft is given "by 

2 . e) 2 0 

• 5 ? ’ 5 ? = 


0 


(2.25) 
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2 

where Cq = G-/ p with G- and ^ as the shear modulus and the 
density of the shaft material respectively, and 0 is the angu- 
lar displacement of the shaft, 

For harmonic motion, the solution of the differential 
equation (2,25) can he written as 

0 (x) = A cos (°Vcq) x + B sin (*">/ c q) x (2.26) 

where fo is the frequency -of oscillation and A and B are the 
constants to he evaluated.. Using the non-dimensional co- 
ordinate "S- , Eqn. 2.26 is written as 


where k 


(2.27) 

1 = non-dimensional torsional frequency -f*- t 


0 = A cos kq$£ + B sin k^ 

CO 

1 - cT 


s' 


The above differential equation can he solved in two parts: 


e. 


Aq cos k qS- + Bq sin kqi 




and 


= 003 ^l^ + B 2 sin k l? ^1^ ^ ^ 1 


(2.28) 


The constants Aq, Bq and AgjBg are obtained from the following 
four boundary conditions: 


(i) 

(ii) 


non-dimensional torque at end A is unity, 

i.e* at 0 = 1 2.29(a) 

torque at end B is zero which gives 

at |=1 4 = ° 2.29(b) 

considering continuity of angular displacement at 


ri 


® L ( fq) - ^ R ( ?q) 


2.29(c) 


(iii) 
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(iv) for the torque balance at the disc 

at ?= = 6 i- 9^ = W A 2 s . © 1/E 2.29(4) 


whero Ijj is the mass moment of inertia of the disc, and I g 
is the mass moment of inertia of the shaft and is equal to 
f I 1 with Ip as the polar moment of inertia of the shaft 
cross-section, 

The (4x4) matrix equation resulting from the use of the 
above boundary conditions is given in Appendix-A. from the 
above set of equations, the constants A’s and B’s are evaluated 
and the receptance is obtained as Q ^ 

Similarly, when unit torque is applied at end B (with 
end A free) the receptances p s -g and p^g are obtained as 6 B 
and 0^ respe ctively* Using the above receptances, the 
propagation constant versus frequency (ih g ) curve is obtained 
using Eqn, 2.9 and the plots are discussed in the next 
chapter,. 


2.3 Evaluation of Receptances for a Periodic Beam Element 

figure 2.5 shows a periodic beam element of length 
1 with a mass Mp at <£• = ^q, where ^ — x/1 with x representing 

the distance of the disc from end A, indicates the rotational 


itiffness of the supports. 

The equation of motion for the beam is given by 


2 <> 4 w 

^ 5 ? 


t) 2 w 

St 5 


= 0 


(2.30) 


where c 2 = with E and I 'as the modulus of elasticity 

and area moment of inertia of the beam,, respectively,, p is 
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the density of team material,. A is the area of cross-section 
of the beam and w is the transverse displacement of the beam, 
lor harmonic motion, the solution of the differential 
equation 2,30 can be written in the non-dimensional form as 


¥(^T ) = cos kg ^ + B 1 s i n C 1 cos - 11 ^2^ 

+ sinh k 2 ^ ( 2- 31 ) 

. 2 

where k^ = 1^ = square of the non-dimensional bending 

frequency, ¥ is the non-dimensional transverse displace-- 

ment = (w/L )/e^ C0 ^, go is the frequency of oscillation with 
A l» C-j_ and as the constants to be evaluated* The 

above differential equation can be solved in two parts s 

¥^ = A R cos k 2 <| + sin k 2 ^ + 0^ cosh k 2 ^ 

+ sinh k 2 1| 0 

¥ R = A 2 cos k 2 ^ + B 2 sin k 2 ^ + 0 2 cosh k 2 ^ (2*52) 

+ D 2 sinh k 2 ^ ^ < r < i 

The constants A' s and B's are obtained from the following 
eight boundary conditions : 


Zero displacement at the ends gives 
(i) at C = 0 \ = 0 


(ii) at 


1 


¥■ 


R 


0 


2.33(a) 

2.33(b) 


Continuity of displacement, slope and bending moment at 
£ = fq gives 
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(iii ) 

at ^ 

= i 

d A = \ 

2.33(c) 

(lv) 

at 

= ? 

H 

II 

nyj 
' — ' 
to 
to 

* 

CM 

(v) 

at 5 

= 5 

_ t ! .1 1 

•1 A = \ 

2.33(e) 

non-di 

.mensional 

bending moment (Ml /El) at end A is 

unity , 

(vi) 

i. e. 

at 

4=0 - s i' + 1 ■ V K = 1 

2.33(f) 



k 1 


at the 

where 

^r 

r 

EI 

is the non-dimensional rotational 

stiffhess /supports 

Moment 

at end B is zero which gives 


(vii) 

at ^ 

= 1 

W R + 75 . £ r . = 0 

2.33(g) 


for the shear force balance at the mass 

(viii) at ^ ^ -%/' + 2 . 33 (h) 

Ki 3 ■■ 

where fl = v .60 is the non-dimensional frequency of 

the beam, is the mass of the beam (=^ll) and M-^ is the mass 
of the disc, 

Ihe ( 8 x 8 ) matrix equation resulted from the above boundary 
conditions is given in Appendix-13. In the some way as described 
in Section 2. 2, the receptances , j3g-g and {3^-g are evaluated from 
the constants A’s and B’s, * Using these receptances, the 
complex propagation constant versus frequency ) curve is 

obtained using Eqn, 2.9 and the plots are discussed in the 
next chapter. 



CHAPTER 3 


RESULTS AND DISCUSSIONS 

3.1 Torsional Wave Propagation Characteristics 

Eirst of all let ns consider an infinite periodic 
shaft shown in Pig. 3«1» The periodic element, IB, of 
this shaft is already shown in Pig. 2.4. Using the analysis 
presented in Chapter 2, the curve of torsional propagation 
constant P- t (= n r t + iM-p p) versus the non-dimensional 
frequency, n g , for this shaft is shown in Pig. 3.2. Both 
the real and imaginary parts of ^ are shown in the same 
figure. The system under consideration Being semi-definite, 
the first propagation zone starts from zero frequency as 
expected. Eurthermore, there are alternate zones of propa- 
gation and attenuation, figure 3.3 shows a similar plot 
for a different value of IpA s with all other parameters 
being the same as in Pig. 3.2. With Ijj equal to zero, i.e. 
for a homogeneous infinite shaft the same curve is as shown 
in Pig. 3.4. As is obvious all frequencies can propagate 
in an infinite homogeneous shaft with no periodic disconti- 
nuity, Equation C-l in Appendix-C explains the nature of 
variation of p ijt withjfn s and is depicted in Pig. 3-4. 

If one consider the other limiting case of 00 * T* e - 

the shaft inertia is neglected, then no propagation zone 
exists as the system becomes an infinite discrete system. 

The propagation constant curve only shows sudden jumps 






FIG- 3-3 VARIATION OF PROPAGATION CONSTANT WITH TORSIONAL FREQUENC 
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in values of from 0 to % and n to 0 whenever p r .j. = 0. 

Hie curves of propagation constant versus frequency 
are indep endent of the variable ^ , Equation G-2 clearly 
shows p^ to depend only on the frequency parameter, i"L s , 
and the inertia ratio Ij)/I g . This result is also to he 
intuitively expected, because the bearings donot play any 
part so far as the torsional oscillations are concerned. 

So long the rotors are at a distance 1 apart, the propaga- 
tion constant curves (Eigs. 3*2 - 3.4) remain unchanged, 

Figure 3.5 shows the band-widths of the first three 
attenuation zones for various values of the inertia ratio 
Ij/tg. Hie band-widths are observed to increase with 
increase in inertia ratio but the rate of increase is more 
for IpAg 41 and it gradually flattens out for higher 
values of I-q/I s . Table 3.1 shows the bounds of the first 
three propagation bands for various values of Ij)A g . The 
propagation band is seen to start from _Q. = (N-l) % . 

As obvious from the discussion presented earlier for the 
two limiting values of 1-^/1 a , Table 3.1 also confirms that 
the propagation band-widths gradually diminish with 
increasing value of 

3.2 Flexural Wave Propagation Characteristics 

We now consider the results of £L exural wave propaga- 
tion constants (p^) for the infinite shaft shown in Tig. 3.6. 
The periodic element of this structure is shown in Tig. 2.5. 
Figures 3.7 - 3.10 show the variation of p^ with frequency 
for various parameters, namely, the mass ratio, -the 
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FIG-3-5 VARIATION OF ATTENUATION BAND-WIDTH WITH INERTIA RATIO 



TABLE 3.1 

BOUNDS OE PROPAGATION BANDS 


1 

v*. 

| Propagation Band 

j First | 

| Second | 

Third 

0.0 

0 0 

N T I N U 

0 U S 

0.1 

0.00 - 2.86 

3.14 - 5.74 

6.28 8.61 


(2.86) X 

(2.60) 

(2.33) 

0.2 

0.00 rr 2.64 

3.14 - 5.32 

6.28 - 8.08 


(2.64) 

(2.18) 

(1.80) 

0.8 

0.00 TT 1.88 

3.14 * 4.22 

6,28 - 6,98 


(1.88) 

(1.08) 

(0.70) 

1.2 

0.00 T 1.62 

3.14 - 3-96 

6.28 - 6.78 


(1.62) 

(0.82) 

(0.50) 

2.0 

0.00 1.32 

3.14 - 3.68 

6, 28 - 6,60 


(1.32) 

(0.54) 

(0.32) 

3.0 

0.00 - 1.10 

3.14 -r 3.52 

6.28 - 6.50 


(1.10) 

(0.38) 

(0,22) 


Quantities in brackets show the widths of the propagation 


hands. 
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disc location parameter, ^ , and the non-dimensional rotational 
stiffness, k r .' for flexural vibration, the system being ; 
positive-definite, the first frequency range 'is always an, 
attenuation zone. 

Comparing figs. 3.7 and 3.8, it is seen that with 
increasing value of ^ (othgr parameters remaining same), both 
the bounds of the odd numbered propagation zones, i.e. first, 
third «... , shift to lower frequencies. Both the bounding 
frequencies of the even numbered propagation zones, however, 
shift to higher values. The difference in this behaviour 
in comparison to the torsional oscillation should be noted. 

This is because the bearings (consequently their locations) 
do affect the manner of flexural vibrations. 

From the figures 3.8 and 3*9, it is observed that 
higher values of the mass^ratio, depress both the 

bounding frequencies of every propagation band. The effect, 
in terms of percentage values, decreases with successive 
propagation tands. • Comparison of lig. 3.7 with Fig, 3.10 
reveals that the rotational stiffness of the bearings 
increases only the lower bounding frequency of every propa- 
gation band where as the upper bounding frequencies remain 
unaffected. This is due to the fact that at the upper bounding 
frequency of every propagation band, the flexural wave mode 

has zero slope at the supports (Ref. 11 ). 

Tables 3.2 and 3.3 show the bounds of the first three 
propagation ’bands for various values of the parameters. It 
is obsarved that for E= 1/2. the mass ratio nanometer does not 
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TABLE 3.2 

BOUNDS OE PROPAGATION BANDS IN ELEXURE 
k r = 0.0 



CQ 

* 1/3 9.9 - 22.4 9.2 - 20.8 7.4 - 16.2 6.2 - 13.2 3.3 - 6.8 

£ ( 12 . 5 ) ( 11 . 6 ) ( 8 . 8 ) ( 7 . 0 ) ( 3 . 5 ) 


1/4 9.9 -' 22.4 

( 12 . 5 ) 


9.4 -- a. 6 

( 12 . 2 ) 


8.0 -' 18.2 

( 10 . 2 ) 


6.9 - 15.1 

( 8 . 2 ) 


3.9 - 8.0 
( 4 . 1 ) 


1/2 


39.5 -" 61.7 39.5 - 61.7 


a 

( 22 . 2 ) 

( 22 . 2 ) 

( 22 . 2 ) 

( 22 . 2 ) 

( 22 . 2 ) 

1 V3 

39.5 -" 61.7 

37.3 -" 56.8 

34.2 -" 48.7 

33.3 -’ 45.5 

32.3 - 41.4 

CO 

( 22 . 2 ) 

0 - 9 . 5 ) 

( 14 . 5 ) 

( 12 . 2 ) 

( 9 . 1 ) 

1/4 

39.5 - 61.7 

36.. 3 - 5 6. '2 

JL . 3 -- 45.9 

29.6 - 41.4 

28c 2 - 35.6 


( 22 . 2 ) 

( 19 . 9 ) 

( 14 . 6 ) 

( 11 . 8 ) 

( 7 . 4 ) 

1/2 

88.8 - 121.0 

82.1 - 112.1 

71.8 - 100.0 

67.9 - 95.7 

63.2 - 90.1 


( 32 . 2 ) 

( 30 . 0 ) 

( 28 . 2 ) 

( 27 . 8 ) 

( 26 . 9 ) 

r c$ 

h 1/3 

88.8 - 121.0 

88.8 - 120.1 

88.8 - 117.9 

88.8 - 116.5 

88 . 8 - 114.7 

Eh 

( 32 . 2 ) 

( 3 !. 3 ) 

( 29 . 1 ) 

( 27 . 7 ) 

( 25 . 9 ) 

1/4 

88.8 - 121.0 

86.0 - 112.9 

83.3 - 102.5 

82.6 -‘ 99.2 

81 . 9 - 95.5 


( 32 . 2 ) 

( 26 . 9 ) 

( 19 . 2 ) 

( 16 . 6 ) 

( 13 . 6 ) 


x Quantitites in Brackets slow tie widths of tie propagation 
lands. 
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TABLE 3.3 


BOUNDS OE PROPAGATION BANDS IN ELEXURE 
k r = 4.0 


g i 

•H j 

-P nd 
03 pjJ 
bn ctf 

] 

hi 

[ 

Disc to 

Shaft Mass Ratio , 


■ 

u 

| 0.0 

0.1 

0.5 

1.0 

5.0 

•P 

CO 

U 

pq 

1/2 

12o8 - 22.4 

(9.6)* 

11.6 - 20.0 
(8.4) 

8.9 - 14.8 
(5.9) 

7.3 - 11.9 
(4.6) 

3.8 - 6 JO 
(2.2) 

1/3 

12.8 - 22.4 
(9.6) 

11.9 - 20.8 
(8.9) 

9.6 - 16.2 
(6. 6) 

8.0 - 13.2- 
(5.2) 

4.3 - 6.0 
(2.5) 


1/4 

12.8 - 22.4 
(9.6) 

12.2 - 21.6 
(9.4) 

10.4 - 18.2 
(7.8) 

9.0 - 15.1 
(6.1) 

5.1 - 8,0 
(2,9) 

■rJ 

cj 

1/2 

42.9 -'61.7 
(18.8) 

42.9 -'61.7 
(18.8) 

42.9 -"61,7 
(18.8) 

42.9 - '61.7 
(18.8) 

42,9 -61.7 
(13.8) 

r~1 

o 

o 

CD 

CO 

1/3 

42.9 - '61.7 
(18.8) 

40.3 -'56.8 
0-6.5) 

36.7 -'48.7 
0-2. °) _ 

35.5 -'45.5 
(10.0) 

34,2 -41.4 
(7*2) 


1/4 

42.9 -"61.7 
(18.8) 

39.3 -"56.2 
(16.9) 

33.5 -'45.9 
(12.4) 

31.4 - 41.4 
(10.0) 

29.5 -35.6 
(6,1) 

nd 

1/2 

92.4 -121.0 
(28.6) 

85.6 -112.1 
( 26.5) 

75.3 -100.0 
(24.7) 

71.4 -"95.7 
(24.3) 

66c 8 -90»1 

(23.3) 

P 

R 

1/3 

92.4 -121.0 
(28.6) _ 

92.4 -120.1 
(27.7) 

92.4 -117.9 
(25.5) 

92.4 -116.5 
(24.1) 

92,4-114.7 

(22.3) 


1/4 

92.4 -121.0 
(28.6) 

89.1 -112.9 
(23.8) 

85.9 -102.5 
(16.6) 

85.0 - ‘99.2 
(14.2) 

84.2 -95.5 
(11.3) 


^Quantities in Brackets indicate the widths of the propagation 
hands « 



41 


affect the bounds of the second hand. This is because at 
these frequencies, the mid-point of every span is a nodal 
point. In fact, the same observation holds good for every 
even numbered propagation bands. For ^ =£ 1/2, all the 
propagation band-widths decrease with increasing mass ratio * 
for the same values of ^ and the propagation band- 

width increases for the successive bands. 

Figures 3*11 and 3*12 show the variation of attenuation 
band-widths with the mass ratio for various values of £ and k r0 
For symmetrical disc location (^ = l/2), the widths of odd 
numbered attenuation bands decrease with increase in mass ratio 
where as these of the even numbered band-widths increase. On 
the , other hand, for unsymmetrical location of the disc, only first 
attenuation band-width is observed to decrease with increase 
in mass ratio and the third attenuation band- width shows an 
increasing trend. In this case, the second attenuation band- 
width shows an initial decrease but later increases at a slow 
rate with increase in mass ratio. However, in all the cases the 
attenuation band-widths do not change significantly beyond Mjj/M-^l* 

3*3 Natural Frequencies of a Finite Shaft 

Figure 3*13 shows an N span finite’ shaft which is 
loaded and supported periodically. The rotational stiffnesses 
at the end supports are taken as half of that at the inter- 
mediate bearings. This is done in order to maintain the 
periodicity of the overall structure. In this section, we 
discuss the the natural frequencies of the shaft both in bending 
and torsional modes. Figure 3*14 shows the variation of 
|i. + versus il . The same figure also contains the plots 

lj u S 




FIG- 




(Qu ut) 

i n niAA-ni inn imimniiTinw 


FIG- 3-12 VARIATION OF ATTENUATION BAND WIDTH WITH MAS 






Torsional frequency parameter,n. s 

F16.3-14 LOCATION OF NATURAL FREQUENCIES(TORSIONAL) 
IN ATTENUATION BANDS 


46 


oj. Abs (l/ft ^ ) in torsion, for various values of ^ . 

As discussed in Chapter 2, in general, ever y propaga- 
tion zone contains ( IT— 1 ) natural frequencies and every 
natural frequency is given "by Abs (l/P^j/ = 0. Thus, out of 
every N natural frequencies, only one frequency depends on the 
value of ^ . These particular natural frequencies depending 
°n ^ are nothing but the natural frequencies of the periodic 
clement shown in Tig. 2.4. Tor symmetric location of the 
disc, i.e., >“ = 1/2, Abs goes to zero at the bounding 

frequencies of the propagation bands. Hence, in this case 
all the natural frequencies can be considered to lie in the 
propagation bands. The natural frequencies lying in the 
propagation bands can be easily determined by the graphical 
procedure explained in Tig. 2.3. These are the values of n 
corresponding to = Hm/U, m varying from 1, 2, . . . . (H-l) . 

The natural frequencies lying in the first three attenuation 
bands are tabulated in Table 3.4. 

Tigure 3.15 and 3.16 show the results for the flexural 
oscillation for two different values of ^ . Again for %= l/2, 
it is seen that all the natural frequencies are contained within 
the propagation zones. Tor ^ 1/2, one out of every U 

natural frequencies lies in the attenuation bands correspon- 
ding to Abs (1 /§ = 0 in bending. Unlike ir. the case of . 
torsional oscillation, in this case the natural frequencies 
lying within propagation bands, also depend on the value of ^ . 
This is due to the dependence of the propagation constant 
itself on the value of § . In every case, the natural 
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TABLE 3.4 

"tii 

N TORSIONAL NATURAL FREQUENCE EOR 
VARIOUS I-j/Eg -AND ^ 


ll 

i VL 


Attenuation bands 


| Eirst 

Seoond 

Third 


0.1 

3.14 

5.74 

9.42 

1/2 

CM 

• 

o 

3.14 

5.32 

9.42 


0.8 

3.14 

4.22 

9.42 


0.1 

3.06 

6.12 

8« 626 

1/3 

0.2 

3.00 

5.96 

8.14 


CO 

• 

o 

2.76 

5.34 

7.36 


0.1 

3.00 

6.27 

8.90 

1/4 

0,2 

2.90 

6.27 

8.40 


0.8 

2.56 

6,27 

7.16 
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frequencies lying within the propagation hands are again 
determined hy the graphical procedure. Table 3.5 gives the 
natural frequencies lying in the first three attenuation hands 
for various values of the parameters involved, 

3,4 Frequency Range of Operation 

We have already seen that mo st of the natural frequen- 
cies of an N span shaft lie within the propagation hands of 

- J* 

the corresponding infinite shaft. In fact, for symmetrical 
location of the disc between the adjacent hearings, all the 
natural frequencies lie within the propagation hands. Hence 
if the torsional and the bending propagation hands are overlapped 
hy proper choice of parameters, a speed range of operation 
can he found which is free from both torsional and £L exural 
resonances. This possibility is explored in this section oy 
attempting to match the propagation hands as far as possible. 

In tliis connection the first torsional propagation hand starting 
from zero frequency is not taken into consideration. The 

bending propagation hand is attempted to he matched with 
the (M+l)^ 31 torsional propagation band. In the discussion to 
follow, symmetrically placed rotors (%= l/2) are assumed to 

be discs of diameter D so that Ijj = / 8 « 

The bending propagation hands are always found to be 
much wider than those in torsion. So as a compromise, instead 
of good matching of the two types of propagation bands, an 
attempt is made to contain the entire torsional propagation band 
within the bending propagation band. It can easily be shown 
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TABLE 3. 5 

N th BENDING NATURAL FREQUENCY EOR 




VARIOUS 

M^/% AND § 




2 r = 

= 4.0 



j M iA 

1 

Attenuation "bands 


5 : 

| First 

Second 

Third 


0.0 

12.79 

42.85 

92.40 


0.1 

11.63 

42.85 

85.55 

1/2 

0.5 

8,925 

42.85 

75.25 


1.0 

7.2466 

42.85 

71.40 


0.0 

12.79 

42.85 

92.40 

1/3 

0.1 

11.913 

40.20 

92.3945 

0.5 

9.560 

35.80 

92.378 


1.0 

7.92 

33.95 

92.370 


0.0 

12.79 

42.85 

92.40 


0.1 

12.209 

39.325 

88.80 

1/4 

0.5 

10,384 

32.95 

84.20 


1.0 

8,89 

30.25 

82.75 
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that for a material with Poisson's ratio = 0.3, /l. ggad m ^ 

D S 

arc related as follows t 

= ^ • rh ,L/a (3 - 1) 

Keeping in mind that the (M+l)^ 11 torsional propagation hand 

always starts at h g = Mm and lies within the range Mm ^ -H. <. 

s 

(M+l) m, if one attempts to match the start of this hand with 
tii 

that of the M bending propagation hand, the resulting values 
of both l/d and D/d are rather low. Much more realistic values 
of these parameters are obtained if one attempts, instead, to 
match the end of these hands. In figures 3.17 and 3.18 both 
the torsional and the bending propagation bands are plotted 
against the common abscissa fl i.e. the bending frequency 
parameter. These two figures refer to two different vaLues 
of the mass ratio parameter M^/M^. i n figures 3.17 ( a ), (b) 
and ( c) two torsional propagation bands are indicated for 
two different sets of values of L/d and D/d. The first set 
of these values matches the start of the M^ bending band 
with that of the (M+l) torsional band, where-as the other 
matches the ends of these bands. The first set of values 
obtained from Pig. 3.17 are also included in Table 3*6. It 
is evident that the values of L/d and D/d are not acceptable 
from practical considerations. In Pig. 3.18 only the ends of 
the bending and the torsional propagation bands are matched. 

With a fixed value of D/d = 3 ■•5, the best possible matching 
obtained in various cases are reported in Table 3.7. In each 
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Mo./Mb= 0-1 

Best matching 

^=l/2,k r =0 ’° 

True matching 

itt'DR 


L/d= 3*51 
D/drl-00 


L /d - 4*63 
D/d=3*50 


L/d=4-79 

D/d=16-7 


Id/Is=27-9 
Extreme case 


Bending frequency parameter, A 5 

(c) ' -■* ■ - 

FIG- 3-1 7 MATCHING OF BENDING AND PROPAGATION BANDS 
(a) M=1 ; (b) M=2 j and (c) M = 3 











jM (Imaginary part) 


a 

c - 2-0 

O) 

a 

£ 


^ -3*0 


L/d=3*92, 

D/d=3-50 


L/d=1*46 

D/d=3*50 


61*6 5 A 
61*70 


15 40 


Bending frequency parameter, n ^ 

' (b) 


Mq /Mb = 1*0 
= 1/2 
*<r =0*0 
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1A3LE 3.6 

PROPAGATION BANDS IN FLEXURE AND TORSION 

V\ = °*1» ?= 1/2, k r = 0.0 


Propaga- 
tion band 

| End frequencies j 
in bending 

| (in-o-fc) | 

End frequencles| L/d > 
in torsion | ratio j 

J (in n b ) | | 

D/d 

jratio 

First 

9.00 - 20.00 

8.997 - 15.816 1.155 

1.22 

Second 

39.45 - 61.70 

39.445 - 54.143 2.532 

loOO 

Third 

82.05 -112.10 

82.046 -100.510 -3.511 

1.00 




TABLE 3.7 


PROPAGATION BANDS IN ELEXURE AND TORSION 


5= 1/2, k r = 0.0 


1 

»jAj 

| Propaga- 1 
tion band j 

| i 

\ J 

| End frequen- 
cies in 
[ flexure 

J (in -*- D ) 

; F~ j 

f Best matching | > 

end frequen- | L/d > 

| cies in,. \ x natio j 

| torsion ( in | ; 

r* 

| D/d 
ratio 

1 


Eirst 

9.00 - 20.00 

15.890 - 19.939 2.04 

3,5 

OJ 

Second 

39.45 - 61,70 

57.174 - 61.544 3.67 

3.5 


Third 

82.05 -112.10 

108.194 -112.100 4,63 

5,5 


Eirst 

6.95 - 14.80 

13.865 - 14.748 1.78 

3,5 

0.5 

Second 

39.45 - 61.70 

60.446 - 61.601 3.88 

3»5 


Third 

71.75 - 99.95 

99,081 - 99.922 4.24 

3.5 


Eirst 

5 » 65 11 • 85 

11.372 - 11,807 1.46 

5=5 

1.0 

Second 

39.45 - 61.70 

61,069 - 61.652 3.92 

j --''D 


Third 

67.90 - 95,70 

95.108 - 95.512 4.07 

3=5 



57 


case, it is seen that the entire torsional hand is contained 
within the corresponding flexural hand. Thus the entire 
attenuation hands for flexural oscillations are free from hoth 
torsional and bending resonances. 

It is already seen in section 3.2, that the start of 
every flexural propagation hand can he shifted to higher 
frequencies without affecting its upper hounding frequency 
hy providing rotational stiffness at the supports. In other 
words, the flexural attenuation hand-widths can he widened 
in this manner. Prom the results presented in I able 3.7, it 
may he noted that there is enough provision for shifting the 
start of the bending propagation hands to higher frequencies 
without interfering with the torsional hands. Hence, upon 
introduction of rotational stiffnesses, k , at the supports, 
a still wider range of operational speed for the shaft can 
he obtained without causing any kind of resonance. 

3.5 Effect of Bearing flexibility in Transverse Direction 

Prom the discussion presented in previous sections, 
it is seen that in the context of wider operating speed 
range of the shaft, it is better to reduce the propagation 
hand-widths in flexural oscillations. To this end, we consider 
an infinite shaft loaded and supported periodically as shown in 
Pig. 3.19. The only difference in this case with that consi- 
dered earlier, is the presence of hearing flexibility, k^. The 
stiffness of the hearing in transverse direction, k^., is 
non-dimensionalised as 
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k t = k t 1 3 /EI (3.2) 

With the introduction of k^, the system "becomes a 
bi-coupled periodic system wherein two pairs of propagation 
constants are involved [ 14 ]. These propagation constants can 
"be determined "by the same receptance method outlined in 
Chapter 2, after modifying the necessary "boundary conditions 
and the equation for the propagation constant. The mathema- 
tical analysis for this purpose is included in Appendix-D. 

The propagation constant curve for a typical case 
is shown in Tig. 3.20. Comparing the first propagation "band 
in this figure with the result given in Table 3.3 for the 
similar case with the following observation can be 

made. 

With the introduction of the bearing flexibility in 
transverse direction, the start of the band remains unaffected, 
whereas the end of the band is shifted to lower frequency. In 
other words, the first propagation band is narrowed down 
considerably. Thus, it may be worthwhile to explore the 
effectiveness of bearing flexibility in widening the operating 
speed range of the shaft. 



Propagation constant 
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FIG.3-20 BENDING PROPAGATION CONSTANTS FOR SHAF.1 
WITH FLEXIBLE SUPPORTS , 





CHAPTER 4 


CONCLUSIONS 

On the basis of the results presented in the thesis, 

several important conclusions can he made. They are listed 

below : 

(i) The wave approach is most suitable for determining 
the bending and the torsional natural frequencies of 
a shaft which is loaded and supported periodically. 
Unlike in other methods, here the computational 
effort remains independent of the number of spans. 

(ii) For an N span periodic shaft, (N-l) torsional natural 
frequencies out of every N remain unaffected by the 
location of the lumped mass within each bay . 

(iii) The attenuation band-widths in torsion do not change 
significantly for inertia ratio greater than unity. 

This is also true in flexure with mass ratio exceeding 
one. 

(iv) By proper choice of parameters within a feasible range, 
(M +l)^ torsional propagation band can be entirely 
contained within the M bending propagation band. 
Perfect matching of these bands yield impractical 
values for the parameters involved, 

(v) The entire flexural attenuation band can be made free 
from both bending and torsional resonances, and hence 
can be chosen as the operating speed range. This 
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range can be effectively increased by providing long 
(introducing rotational stiffness k r ), transversely 
flexible (introducing k t ) bearings. 
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APPENDIX - A 


MATRIX EQUATION EOR THE BOUNDARY CONDITIONS 
GIVEN BY EQUATIONS 2.29 (a-d) 


0 

1 

0 

0 

t \ 

A 1 


'l/K 

0 

0 

-kq sin k^ 

kq cos 

h 


0 ’ 

COS 

sin kq 

-COS kq ^q 

-sin kq )^q 

r 2 

f 

° > 

~kq sin kq ^ q 
^D^s* 

k-^ cos ^ 

- j bA s * a I* 

kq sin kq ^ q 

-kqCoskq<jq 

i - 
1 

ki 
! : 

: 0 

COS ^ *j^ 

sin k-^^q 


- 

\ * 

1 

(A.l) 



Erom the above matrix equation, the constants Aq, Bq, A 2 and 
Bg can be easily determined. 

Similarly, when the non-dimensional torque at end B is unity, 
boundary conditions 2.29 (a) and (b) change to 

at ^ = 0 e’i = 0 

and at ^ = 1 O f R = 1 

Thus the right hand matrix in Eqn. A.l gets modified to 

( \ 

0 

1 


0 



Prom the above matrix equation, 

A 2 1:0 ?2 are determined and the 
a s W^, 


the constants A 1 to and 
receptance p^A is obtained 


Similarly, when unit non-dimensional moment (Ml /El ) is 
applied at end B f the boundary conditions 2.33 (f) and (g) 
change to i 



at $ 

= 0 

¥g + - 0 

and 

at $ 

= 1 

— t » q _ 

Wg + r' w R = respectively. 

Then the 

right 

hand matrix 

in Eqn. B,1 changes to 


0 

0 

0 

0 



0 

-1 

0 


A new set of values for the constants Aq-D-^ and A 2 *-D 2 are 
determined using the above right hand side matrix. Then the 
receptances p-g-g and (3^ are obtained as ¥g and respectively. 


APPENDIX - 0 


PROPAGATION CONSTANT IN TORSION 

C-l Propagation Constant for a Homogeneous Shaft 

Prom the matrix equation of Appendix - A the constant 
A's and B’s are evaluated for Ij)/l g . = 0,0 and the receptances 
§AA ^BA are obtained as 


Paa - e 


h + P 2 

U 3 + 


! fc, 3c, 3c. 

where u i - - ^ sin ^-(3^ sin \ tan cos ~^) 


I D - 1l 1 lz l ^ 

Ug = + (j— . k-^.sin Y" ~ cos ^-)(cos -g=-+ sin tan k B ) 


kp k^ k-, 

U 3 = cos 7j~ (k-j_ sin ^ lc B tan k-j_ cos -*r) 

k*L 1 B 2 k. k-| 3c. 

= — (k-^ sin cos Ty— ) (cos + sin tan k B ) 


and p BA = 6 B = [l/cos kj- 1/(U 5 + U^)] 


P 


Ail 


Hence, cosftii^. = *p — - = cos k^ = cos JX 


or. 


BA 

^i,t = ( A e( 


( 0 . 1 ) 


C-2 Dependence of Propagation Constant on 

Prom the first set of boundary conditions (2.29(a~d)) , 
where the torque is applied at end A, receptances and p B ^ 


are found to be 



C-2 


&AA “ ® A 

and p BA = & b 

Similarly from the second set of boundary conditions, 
when -unit torque is applied at end B receptances and 
are obtained from appendix A, as 

■ U g 

^BB = ®B = tig 

ry 

where Ug = Qos^k-^^ .sec k^sin fc^ 

+ sin k-g tan k-j_ tan k^j 2 - tan kj|- sec k^)] - 1 

U 6 = k x tan k x + X-g/l s .k|. cos 2 k^ (l+tan k^ tan k x ) 

o o 

and = - l/[ co s k-j_ | k^ tan k-j_ + I-g/lg.k^.cos k^ (1 

+ tan k^ tan 

Substituting the above four receptance functions in Bqn. 2.10, 
it is .seen fhat the torsional propagation constant, can 

be obtained from 

, ' P BB ^AA T A B1 - n A 

cosh p t = — — = cos St s - Ij/ig- ’y* sinu s 

H AB 

So, the propagation constant, p^, is independent of 


-L *r\ /y 

= [l - { j—.k^.sin k^ (cot k-^lj + tan k^)}]/ 

3 

I ys ry ry 

[k^ tan k^ + j"*k^. cos (l+tan' k^ tan k^)] 
s 

I ■ 

= l/[cos k^ { k-j. ba 21 ic i + j^.b; 2 . cos 2 k 1 ^(l 
+tan k-j^ tan 1^)^] 



APPENDIX - D 


EQUATIONS WITH TRANSVERSELY ELEXIBDE BEARINGS 

As now there are two coupling co-ordinates at either end 
of every elejnent, the system under consideration is a hi-coupled 


system. 

Shear 

force 

balance 

at 

both the ends gives 


(i) 

at ^ = 

0 

\ A t 

+ 

iy>|H 

II 

o 

(D.l) 

(ii) 

at T f = 

1 

- W R At 


I*r = 0 

(D.2) 


Continuity of displacement, slope and bending moment 
at gives 

(iii) at ^ = ^ ] _ W L = W R (D.3) 

(iv) at |^= 5 1 Wj? = W R (D.4) 

(v) at ^ ¥^’ = W R (D.5) 

Non-dimensional bending moment (Ml/El) at end A 
is unity 

(vi) i.e. at ^ = 0 -Wq + = 1 (D*6) 

Moment at end B is zero which gives 

(vii) at ^ = 1 i" + = 0 (D.7) 

Por the shear force balance at the mass 

( viii ) at ^ (D.8) 


The (8x8) matrix equation resulting from the above 
boundary condition is as follows : 
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D -3 

In this analysis, the receptance equations given in Chapter 2 
for a mono-coupled system, are modified as : 


11 

<b 

P AA 

m a 

+ 

P AB 

“b 

+ 


J A 

4- 

^AB 



^ - 
B 

P BA 

m a 

+ 

^BB 

“b 

+ 

^BA 


4- 

^BB 

*B 

(D. 10 ) 

II 

<r| 

1^ 

6 AA 

m a 

+ 

6 AB 

m b 

~h 

a AA 

t a 

4* 

a AB 

P B 


II ' 

6 ba 

m a 

+ 

6 bb 

m b 

+ 

a BA 

P A 

+ 

a BB 

*B 



where is the harmonic rotation at i to 'unit harmonic 
moment at j, 

j is the harmonic rotation at i due to unit harmonic force at j, 
6^^ is the harmonic deflection at i due to unit harmonic moment at j, 
and a^-j is the harmonic deflection at i due to unit harmonic 
force at j. 

four sets of receptance functions, namely, 

6 AA and S BA ], [p AB , P BB , 6 iiB said 6 BB ] s [ i BA , oc u and 

and [ ^3, >33 , a^-g and CC33] are obtained respectively 
by four sets of boundary conditions which are resulted by 
the application of unit moment and unit force, respectively, 
at each end. The resulting four right hand column matrices 
are as follows: 
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Unit moment Unit moment Unit shear force Unit shear force 
about A about B at A at B 


0 

0 

0 

0 

0 

1 
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0 

0 
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